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SECTION  1 
INTRODUCTION 


BACKGROUND 

The  present  AFWL/DES  ground-shock  and  cratering  computer  code  (AFTON)  has  in¬ 
ternal  difficulties  which  limit  its  ability  to  perform  calculations  of  large, 
turbulent  deformations  in  regions  in  which  material  mixing  and  material 
strengths  are  significant.  This  limitation  has  led  to  severe  problems  when 
attempting  to  model  late- time,  earth  material  response  to  high-explosive 
detonations.  This  problan  can  be  remedied  by  modifying  the  CERF  version  of 
the  Large-Displacement,  Elastic/Plastic  Computer  Code  (LDEPCC).  This  code, 
CSQII,  originally  developed  to  calculate  large  deformations  resulting  from 
underground  nuclear  detonations,  has  been  applied  to  cratering  by  hypervel¬ 
ocity  impacts  on  steel  and  aluminum  plates. 


OBJECTIVE 

The  purpose  of  this  effort  was  to  modify  the  CERF  version  of  LDEPCC  so  as  to 
make  it  compatible  with  the  AFWL/DES  ground-shock  and  cratering  calculational 
needs,  and  thereby  enable  the  Air  Force  to  Increase  its  knowledge  and  under¬ 
standing  of  nuclear  effects  on  hardened,  protective  structures. 


SCOPE 

After  the  performance  of  a  test  case  to  check  the  validity  of  the  results  of 
the  code  against  an  available  exact  solution,  the  material  models  derived  from 
Cylindrical  In-Situ  Test  (CIST)  15  were  to  be  used  to  compare  the  calculated 
and  observed  large-deformation,  earth  material  response  of  the  Pre-Dice  Throw 
100-Ton  TNT  Event.  This  evaluation  could  not  be  performed  because  of  diffi¬ 
culties  in  the  treatment  of  mixed  material  cells;  however,  a  special  test  was 
performed  as  a  partial  evaluation. 
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SECTION  2 

CSQ  II  COMPUTER  CODE 


CSQ  II  is  a  Fortran  program,  written  by  Dr.  Samuel  L.  Thompson  of  Sandia  Labor¬ 
atories,  for  the  computation  of  two-dimensional  material  motion  in  either  a 
rectangular  or  a  cylindrical  coordinate  system  (ref.  1).  A  Eulerian  finite- 
difference  method  is  employed  with  a  high  degree  of  automation  built  into  the 
program  string.  Simultaneous  energy  flow  resulting  from  radiation  and  thermal 
conduction  can  be  treated  as  well  as  a  variety  of  solid  dynamic  properties. 

Many  aspects  of  CSQ  II  are  similar  to  those  of  the  one-dimensional  code  CHART 
D;  in  particular  the  CHART  D  equation-of -state  is  used  in  CSQ  II. 

Essentially  CSQ  II  is  a  hybrid  of  Lagrangian  and  Eulerian  methods.  During  each 
cycle  of  computation  the  mesh  is  treated  as  materially  fixed,  and  a  complete 
Lagrangian  calculation  is  performed  as  the  mesh  moves  with  the  material.  After 
this  phase  the  code  enters  a  rezoning,  or  Eulerian  section,  in  which  the  entire 
mesh  is  redefined  and  moved  back  to  the  position  it  had  before  the  Lagrangian 
phase. 

A  unique  feature  of  CSQ  II  is  the  retention  of  cell  void  volumes.  This  is  ne¬ 
cessary  for  the  computation  of  solid  material  responses  in  the  Eulerian  method, 
since  edges  of  solids  and  interior  spalling  can  be  properly  handled  only  in 
this  manner. 

The  standard  version  of  CSQ  II,  which  can  treat  up  to  400  meshes  in  either 
spatial  coordinate  system  (160,000  cells),  requires  130,000.  to  270,000  words 
of  central  memory  (CM)  on  the  CDC  7600,  depending  on  the  number  of  different 
materials  in  the  problem  and  the  number  of  special  routines  required  (e.g., 
high  explosives,  energy  flow,  etc.). 


1.  Thompson,  S.  L.,  CSQ  —  A  Tido-Dimenaional  Hydrodyncania  Program  with  Energy 
Flow  and  Material  Strength,  SAND  74-0122,  Sandia  Laboratories,  Albuquerque, 
New  Mexico,  August  1975. 


The  program  retains  the  complete  Identity  of  each  material  (up  to  ten)  present 
In  the  calculation.  For  each  material  over  two,  three  additional  cell  storage 
variables  are  required.  From  9  to  40  variables  per  cell  are  retained,  depend¬ 
ing  on  the  problem;  a  special  preprocessor  program  Is  used  to  tailor  the  stor¬ 
age  to  the  problem. 


The  main  cell  storage  arrays  are  packed  to  avoid  wasted  space.  The  packed 
arrays  are  kept  In  extended  core  storage  (ECS),  with  only  a  small  portion  of 
the  entire  array  held  In  CM  at  any  time.  In  practice,  the  largest  problem 
(in  terms  of  the  number  of  cells)  is  limited  by  the  amount  of  ECS  available. 

Four  optional  subroutine  sets  may  be  used  with  CSQ  II.  The  following  addi¬ 
tional  CM  capacity  Is  required  when  these  optional  sets  are  used: 


Energy  Flow 

Elastic/Plastic 

Energy  Sources/ 

High  Explosives 

Internal  Fixed  Mesh 


350O3 

3400. 

2300g 

5OO3 


The  CSQ  II  package  contains  seven  programs  and  their  general  functions  are  as 


follows: 

(1) 

PRECSq 

Preprocessor.  Tailors  storage  arrays  of  remaining 

(2) 

CSQGEN 

programs.  (Must  be  used  before  any  of  the  remain¬ 
ing  programs.) 

Generator.  Generates  new  problems  and  rezones  re¬ 

(3) 

CSQ 

sults  of  previous  calculations. 

Main  Program.  Runs  problems. 

(4) 

CSQLINE 

One-Dimensional  Plot  Program 

(5) 

CSQPLT 

Two-Dimensional  Plot  Program 

(6) 

CSQSURF 

Tnree-Dimensional  (Surface)  Plot  Program 

(7) 

CSQTa,P 

Output  Data  Tape  Editor 
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SECTION  3 
IMPLEMENTATION 


The  only  modification  required  to  implement  the  basic  CSQ  II  on  the  AFWL  CDC 
7600  System  was  the  addition  of  a  subroutine.  This  subroutine  (HOROLOG),  which 
supplies  date,  time,  and  time  remaining  in  the  computer  run,  is  a  library  sub¬ 
routine  in  the  Sandia  Laboratories  CDC  7600  System  but  does  not  exist  on  the 
AFWL  System.  However,  HOROLOG  was  included  simply  by  combining  AFWL  System 
Library  Subroutines  XTIME,  DATE,  and  TIME. 

Validation  of  CSQ  II  consisted  of  examining  the  velocities  and  displacements 
produced  by  applying  a  constant  pressure  over  a  spherical  surface  buried  in  a 
multilayered  geology.  These  responses  were  compared  to  an  exact  analytic  elas¬ 
tic  solution. 


INTERPOLATION 

During  the  running  of  the  validation  problem,  numerous  large  irregular  discrep¬ 
ancies  were  noted  when  the  CSQ  II  results  were  compared  with  the  analytic  solu¬ 
tion.  These  discrepancies  were  traced  to  the  technique  used  to  determine  tracer 
particle  velocities  within  a  cell.  This  technique  determines  tracer  particle 
velocities  by  interpolating  on  the  velocities  associated  with  the  cell  contain¬ 
ing  the  tracer  particle,  independent  of  the  velocities  in  neighboring  cells. 

This  results  in  errors  when  the  tracer  particle  is  very  near  one  of  the  bound¬ 
aries  of  the  cell  containing  it.  Thus,  a  new  interpolation  scheme  was  written 
to  correct  this  situation.  Details  of  this  scheme  are  presented  in  appendix  A. 
With  this  new  tracer  particle  velocity  interpolation  scheme,  the  agreement  be¬ 
tween  CSQ  II  responses  and  the  analytic  solution  improved  from  clearly  unaccept¬ 
able  to  very  good. 


SYMBOLOGY 

A  minor  change  was  made  to  the  symbology  used  by  CSQ  II  in  printing  out  the 
material  plot  during  the  long  edit  cycle.  The  ten  materials  of  CSQ  II  are 
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designated  to  be  a  single  material,  called  DOT^  and  the  remaining  nine  materials 
are  called  x.  Formerly,  In  the  material  plot  the  DOT  material  was  printed  out 
as  a  single  .  and  the  nine  X  materials  were  printed  out  as  the  digits  2  through 
9  and  the  single  letter  A.  This  has  been  changed  to  print  out  the  nine  X  mate¬ 
rials  as  the  digits  1  through  5,  omitting  the  letter  A.  It  Is  felt  that  this 
Is  more  In  keeping  with  the  mental  picture  one  forms  of  the  ten  materials,  the 
nine  X  materials  in  particular. 


AFWL  ENGINEERING  (AFWLE)  EQUATION-OF-STATE 

The  only  real  problems  encountered  In  the  Implementation  of  CSQ  II  are  related 
to  the  inclusion  of  the  AFWLE  equation-of-state  In  CSQ  II.  Specifically,  four 
probl'^s  were  encountered  in  this  area. 

First,  the  processing  of  the  AFWLE  input  variables  was  greatly  complicated  by 
the  fact  that  the  many  automated  features  of  CSQ  II  require  the  special  setting 
of  numerous  flags  and  variables  at  Input  time.  Relatively  few  of  these  flags 
and  variables  are  documented  as  to  purpose,  value,  or  interaction  with  one  an¬ 
other.  These  had  to  be  discovered  by  tracing  the  path  of  one  of  the  existing 
CSQ  II  equations-of-state.  Unfortunately,  the  eight  equatlons-of-state  built 
into  CSQ  II  do  not  exist  as  eight  distinct  modules,  but  they  are  mixed  together 
In  a  rather  complex  involuted  way  with  numerous  additional  flags  and  pointers 
to  point  the  way  through  this  labyrinth.  The  result  of  all  this  Is  that  there 
is  always  a  nagging  doubt  as  to  whether  or  not  everything  has  been  properly  set. 

Second,  CSQ  II  employs  a  temperature-dependent  equation-of-state,  while  that  of 
AFWLE  is  energy  dependent.  That  is,  whenever  CSQ  II  calls  the  equation-of-state, 
the  calling  parameters  are  density  and  temperature;  AFWLE  expects  density  and 
internal  energy.  The  basic  form  of  the  equation-of-state  used  in  AFWLE  is 

P  =  fj(p)  +  Tpe  (1) 


where 

P  =  pressure 

fi(p)  =  a  given  function  of  the  density 

r  =  the  Griinelsen  ratio  (a  material  constant) 
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p  *  density 
e  *  internal  energy 

The  problem  is  to  determine  e,  given  the  temperature. 


In  CSQ  II  the  following  relationship  between  energy  and  temperature  exists; 

e  =  C^(t  -  t„)+  (2) 


where 

Cy  *  heat  capacity  of  the  material 
t  =  temperature 

t^  =  temperature  on  the  Hugoniot 

=  mechanical  energy  due  to  compression  (or  expansion)  on  the  Hugoniot 
If  is  the  pressure  along  the  Hugoniot,  is  related  to  by 


where 


PhH 


(3) 


n  »  1  -  Pq/p 

The  AFWLE  equation-of-state  is  of  the  Mie-Griineisen  form,  and  can  be  written 


P  *  P^  +  rp(e  -  e^) 

=  P^  -  rpE„  +  rpe  (4) 


From  eq.  (4)  and  eq.  (1)  we  see  that 


or 


f,(p)  =  Ph  -  rpE^^ 


P^  “  fj(p)  +  rpE^ 


Combining  eq.  (3)  and  eq.  (5)  yields 


f,(p)n 
Eu  = 


H  2p^  -  pFn 


In  general,  t^  is  given  by 


Jn 


f  -  f  3.  e  r  «-rri  2  d  . 

-  V  j  e  n 


(5) 


(6) 
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Unfortunately,  the  form  for  yields  an  expression  In  the  Integral  which  can¬ 
not  be  evaluated  simply.  The  approach  employed  here  was  to  use  an  approxima¬ 
tion  which  has  the  correct  limiting  values  and  derivatives  near  p  *  p^. 

where  t^  Is  the  ambient  temperature. 

The  values  for  and  t^,  given  by  eqs.  (6)  and  (7),  respectively,  are  used  In 
eq.  (2)  to  determine  e  for  the  AFWLE  equation-of-state  [eq.  (1)]. 

To  compute  the  pressure  at  low  densities  (p  <  p^)  AFWLE  uses  the  form 

P  =  p[h  +  (r  -  H)  ^  -  e^(l  -  exp  Nn(l  -  n))]  (8) 

where 

U  =  P/Pq  -  1 

Eg  =  sublimation  energy  for  the  materiel 
N  =  Cj/re^ 

Cq  »  bulk  sound  speed 
H  *  Y  -  1  ■ 

Y  »  ratio  of  specific  heats  of  the  perfect  gas 
Since  neither  H  nor  y  Is  specifically  stated  in  the  AFWLE  input,  H  is  set  to 
0.4,  a  value  which  approximates  the  value  for  most  gases. 

In  this  case  (p  <  p^),  e  is  computed  from 

e  =  Cy  (t  -  t,^) 

When  p  =  p^,  eq.  (8)  reduces  to 

p  =  rpe  =  rpC^(t  -  t^) 

and  is  continuous  with  eq.  (1)  at  this  density. 

Third,  when  CSQ  II  determines  the  pressure  within  a  mixed  material  cell,  it  first 
determines  the  pressure  for  each  material  in  the  cell  and  does  a  volume  averaging 
to  obtain  the  net  pressure  for  the  entire  cell.  Having  the  cell  pressure,  CSQ 
II  then  adjusts  the  density  of  each  material  in  the  cell  so  as  to  yield  the 


computed  cell  pressure  for  each  Individual  material.  This  process  of  individual 
material  density  adjustment  is  essentially  a  process  of  Newtonian  iteration  and 
requires  that  the  derivative  of  pressure  with  respect  to  density  be  known  for 
each  material  in  the  cell.  This  derivative  must  also  be  continuous  and  well 
behaved  in  the  region  where  this  calculation  is  to  take  place. 

Unfortunately,  the  form  for  pressure  in  AFWLE  does  not  have  derivatives  which 
are  well  behaved.  In  particular,  the  function  fj(p)  of  eq.  (1)  is  a  piecewise 
linear  function  and  thus  the  derivative  of  pressure  with  respect  to  density  is 
a  step  function.  Similarly,  the  derivative  of  pressure  given  by  eq.  (7)  behaves 
very  erratically  over  much  of  the  range  p  <  p^.  The  result  of  this  has  been 
convergence  failure  when  computing  densities  in  mixed  cells.  Such  convergence 
failure  causes  automatic  run  termination  of  CSQ  II.  Much  effort  has  been  spent 
attempting  to  correct  this  problem  but  without  success.  Indeed,  there  may  not 
be  a  solution  strictly  within  the  framework  of  the  AFWLE  equation-of-state. 

Fourth,  the  AFWLE  equation-of-state  is  history  dependent.  That  is,  the  exact 
form  of  the  equation-of-state  depends  on  the  previous  states  of  the  material 
under  consideration.  The  equation-of-state  behaves  very  much  like  a  porous 
material.  Or  foam.  On  loading  the  material  compresses  elasticly  up  to  a  cri¬ 
tical  pressure,  from  which  point  it  compresses  plasticly.  On  release  it  does 
not  unload  along  the  load  path  but  rather  unloads  along  an  elastic  path  from 
the  point  of  maximum  compression.  On  reloading  the  material  compresses  elas¬ 
ticly  up  to  the  previous  point  of  maximum  compression  where  it  again  continues 
its  plastic  compression.  This  requires  that  the  state  of  maximum  compression 
be  retained  among  the  cell  storage  variables  for  each  material.  Because  this 
would  require  considerable  additional  storage  and  a  complete  rearrangement  of 
the  CSQ  II  storage  scheme,  it  was  decided  that  this  feature  would  not  be  in¬ 
cluded  at  this  time.  The  scheme  implemented  in  CSQ  II  is  to  allow  the  material 
to  unload  along  the  original  load  path.  On  recompression  the  material  loads 
along  the  same  path  as  the  original  compression. 

The  modified  input  required  for  the  AFWLE  equation-of-state  is  presented  in 
appendix  B. 
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SECTION  4 
EVALUATION 


Originally  it  had  been  intended  to  evaluate  CSQ  II  by  comparing  the  results  of 
the  Pre-Dice  Throw  Experiment  with  the  results  of  a  CSQ  II  simulation  of  that 
experiment.  Briefly  stated,  the  Pre-Dice  Throw  Experiment  consisted  of  the 
detonation  of  a  100-ton  charge  of  TNT  over  a  multilayered  earth  medium  and  the 
examination  of  the  velocities  and  displacements  produced  in  the  underlying  me¬ 
dium  at  various  test  points  out  to  a  50-msec  problem  time.  However,  because 
of  difficulties  in  the  treatment  of  mixed  material  cells  (section  3),  no  run 
was  completed  beyond  1.9  msec;  this  is  just  a  little  beyond  the  time  of  com¬ 
plete  detonation  of  the  TNT  and  no  significant  motion  was  imparted  to  the  un¬ 
derlying  earth  by  this  time.  Thus,  no  significant  evaluation  of  CSQ  II  was 
possible  from  these  runs. 

A  partial  evaluation  of  CSQ  II  without  the  AFWLE  equation-of-state  was  made 
by  comparing  the  exact  analytic  solution  of  a  special  test  case  with  the  CSQ 
II  solution  for  that  problem.  The  test  case  consisted  of  a  constant  2000-psi 
pressure  applied  uniformly  over  the  interior  surface  of  a  spherical  cavity  in 
a  layered  geology  (fig.  1).  Of  particular  interest  were  the  shock  arrival 
times  and  the  velocities  and  displacements  produced  by  the  shock  waves  at  vari¬ 
ous  target  points  near  the  interface  between  the  two  materials.  In  general, 
shocks  near  the  interface  are  a  result  of  direct  waves  and  waves  reflected  or 
refracted  at  the  interface.  The  results  of  the  CSQ  II  run,  together  with  the 
analytic  solutions,  for  selected  peri -interfacial  target  points  are  given  in 
figures  2  through  7. 


The  only  problem  encountered  during  the  running  of  the  test  case  was  the  large 
discrepancies  in  velocities  and  displacements  caused  by  the  CSQ  II  interpola¬ 
tion  schemes.  This  problem  is  thoroughly  discussed  in  section  3;  the  new  cor¬ 
rective  interpolation  technique  is  presented  in  appendix  A. 


11 


. 


Analytic  Solution 
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V 


Displacement,  in  Velocity,  in/sec 


Displacement,  in  Velocity,  in/sec 


M 


6.00E+01 


4.00E+01U 


2.00E+01^ 


•2.00E+01 

7.50E-01 


5.00E-01 


2.50E-0 


•2.50E-01 

l.C 


Analytic  Solution 
CSQ  II  Solution 

V-Vertical 

H-Horizontal 


SECTION  5 

CONCLUSIONS  AND  RECOMMENDATIONS 


There  are  two  avenues  of  approach  open  for  the  future  use  of  CSQ  II.  First, 

CSQ  II  could  be  used  in  its  present  form,  i.e.,  without  the  AFWLE  equation-of- 
state.  However,  when  problems  to  be  run  with  CSQ  II  are  defined  in  terms  of 
the  AFWLE  equation-of-state,  as  will  often  be  the  case,  it  will  be  necessary 
to  convert  the  AFWLE  input  parameters  to  a  form  compatible  with  that  built  in¬ 
to  the  CSQ  II  equation-of-state.  This  means  approximating  the  jointed  linear 
AFWLE  equation-of-state  with  a  polynomial.  A  computer  program  could  be  devised 
to  accomplish  this.  The  second  approach  is  to  spend  the  additional  time  and 
money  requi»*ed  to  fully  implement  the  AFWLE  equation-of-state  into  CSQ  II.  Be¬ 
cause  of  the  highly  involuted  logic  flow  of  CSQ  II  and  the  complex  iteration  of 
its  many  automated  features,  this  would  be  no  simple  task.  Communications  with 
Dr.  S.  L.  Thompson  of  Sandia  Laboratories,  the  author  of  CSQ  II,  have  confirmed 
this  assessment.  However,  should  this  approach  be  chosen,  it  would  be  advis¬ 
able  to  contact  Dr.  Thompson  since  he  has  as  much  knowledge  as  anyone  in  the 
internal  complexities  of  CSQ  II. 
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APPENDIX  A 

IMPROVED  TECHNIQUE  FOR  COMPUTING  INTERPOLATED 
VELOCITIES  OF  TRACER  PARTICLES 

Formerly  when  a  tracer  particle,  T  (fig.  Al),  was  located  in  cell  (i.j),  the 
x-component  of  its  velocity,  Vj^(T),  was  calculated  from 

Vx(T)  =  x[v^(i  +  l,j)  +  (1  - 

where 

X  » 

Vjj(i  +  l,j)  *  velocity  (x-component)  of  the  right  boundary  of  cell  (i,j) 
Vj^(i.j)  ■  velocity  of  the  left  boundary  of  cell  (i,j) 

X(T)  =  x-coordinate  of  tracer  particle 
X(i)  “  x-coordinate  of  left  boundary 
X(i  +  1)  =  x-coordinate  of  right  boundary 
The  y-component  of  the  tracer  particle  velocity,  Vy(T),  was  calculated  in  a 
similar  manner. 


Figure  Al.  Tracer  Particle  Velocity  Interpol  cation  in  CSQ  II 


21 


This  technique  led  to  errors  in  the  calculated  velocities  when  the  tracer  par¬ 
ticle  was  near  one  of  the  cell  boundaries.  For  example,  when  the  particle  is 
on  or  very  near  the  lower  boundary  of  cell  (i,j),  it  is  expected  that  the  cal¬ 
culated  particle  velocity  would  be  influenced  as  much  by  the  velocities  of  the 
cell  below,  cell  (i,j  -  1),  as  by  the  velocities  of  cell  (i,j)  itself.  In  the 
former  technique,  however,  the  particle  velocity  was  not  affected  by  cell 
(i,j  -  1)  at  all;  this  led  to  very  irregular  velocity  behavior  when  the  posi¬ 
tion  of  the  particle  fluctuated  from  one  side  of  the  boundary  to  the  other. 

In  the  improved  interpolation  scheme,  when  a  tracer  particle  is  located  in 
cell  (i,j)  its  velocity  is  calculated  from  the  velocities  associated  with 
cell  (i.j)  and  the  velocities  of  one  or  more  of  the  neighboring  cells.  The 
example  presented  in  figure  A2  illustrates  this  new  technique. 

Since  T  is  located  in  the  lower  half  of  cell  (i,j),  the  x-component  of  the  ve¬ 
locity  of  T,  V^(T),  will  be  interpolated  from  the  velocities  of  the  cell  below 
cell  (i,j);  i.e.,  cell  (i,j  -  1).  If  T  is  in  the  upper  half  of  cell  (i,j),  in¬ 
terpolation  would  be  carried  out  with  cell  (i,j  +  1),  the  cell  above  cell  (i,j). 
The  interpolation  is  bilinear;  i.e.,  linear  in  both  the  x-  and  y-directions. 

First  velocities  are  interpolated  in  the  x-direction  in  a  manner  identical  to 
that  in  the  former  interpolation  scheme.  This  yields  for  cell  (i,j)  and  cell 
(i.j  -  1).  respectively, 

V^(T,j)  *  X  V^(i  +  l,j)  +  {1  -  X)  V^(i,j) 

V^(T,j  -  1)  =  X  V^(i  +  l.j  -  1)  +  (1  -  X)  V^(i,j  -  1) 

where 

X  * 

Vj^(i,j)  =  velocity  of  left  boundary  of  cell  (i,j) 

Vj^(i  +  l,j)  =  velocity  of  right  boundary  of  cell  (i,j) 

Vj^{i,j  -  1)  =  velocity  of  left  boundary  of  cell  (i,j  -  1) 

Vj^(i  +  l,j  -  1)  =  velocity  of  right  boundary  of  cell  (i.j  -  1) 

Next,  the  values  V  (T,j)  and  V  (T,j  -  1)  are  interpolated  (linearly)  in  the 

A  A 

y-direction  to  obtain  V^(T). 


> 
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Figure  A2.  Improved  Tracer  Particle  Velocity  Interpolation  for  CSQ  II 

If  Yj  denotes  the  y-coordinate  of  the  center  of  cell  (i,j)  and  ^  denotes 
the  y-coordinate  of  the  center  of  cell  {i,j  -  1), 

=  (^j  *  *  i)/^ 

and 

See  figure  A2.  If 

y(T)  -  V, 

\’= - - 

-  1 

where  y(T)  is  the  y-coordinate  of  the  tracer  particle, 

Vx(T)  »  ^y[Vx(T.j)  +  (1  -  X)  Vx(T,j  -  1^ 

It  should  be  noted  that  this  bilinear  interpolation  is  commutative.  That  is, 
the  results  would  be  exactly  the  same  had  the  first  interpolation  been  in  the 
y-direction  followed  by  interpolation  in  the  x-direction.  A  similar  bilinear 
interpolation  process  is  used  to  determine  Vy(T),  the  y-component  of  the  tracer 
particle  velocity. 
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APPENDIX  B 

MODIFIED  INPUT  INSTRUCTIONS 


Complete  input  instructions  for  CSQ  II  are  given  in  references  2  and  3.  The 
following  modifications  do  not  change  the  basic  CSQ  II  input  but  rather  supple 
ment  that  input  for  the  AFWLE  equation-of -state. 


CARD  SET  3:  ANALYTIC  EQUATION-OF-STATE  DATA 

Card  set  3  (ref.  3)  provides  the  equation-of -state  data  for  the  various  mate¬ 
rials  in  the  problem  to  be  run  with  CSQGEN.  When  the  AFWLE  equation-of-state 
is  to  be  used,  the  following  cards  must  be  included,  one  set  for  each  AFWLE 
type  material.  There  is  no  problem  in  mixing  AFWLE  type  materials  with  regu¬ 
lar  CSQ  II  equation-of-state  materials  in  a  given  problem.  Equation-of-state 
parameters  for  CSQ  II  type  materials  are  given  in  reference  2.  This  card  set 
is  in  the  CSQGEN  input  stream. 

Card  3.1  -  FORMAT  (13,  15,  12,  5A10,  2E10.3) 


Variable 

Columns 

Description 

1 

1-3 

Equation-of-State  Number.  Always  -1  to  -20. 

2 

4-8 

Must  be  zero  for  AFWLE  equation-of-state  read  in  from 
cards;  nonzero  for  library  equations-of-state. 

3 

9-10 

Used  only  by  library  equation-of-state  routines;  set 

to  blank. 

4-8 

11-60 

Identification  Label.  Any  BCD  information. 

2.  Thompson,  S.  L.,  Improvements  in  the  CHART  D  Energy  Flow-Hydrodynamia 
Code  V:  2972/1973  Modifications ^  SLA-73-0477,  Sandia  Laboratories, 
Albuquerque,  New  Mexico,  October  1973. 

3.  Thompson,  S.  L. ,  Input  Instructions  for  CSQ  II  -  A  Ten  Material  Version 

of  CSQ,  RS  5166/135,  Internal  Memorandum,  Sandia  Laboratories,  Albuquerque, 
New  Mexico,  October  1975. 


Variable 

Columns 

Description 

9 

61-70 

RHUG  -  Initial  Density  for  Hugoniot  Calculation.  If  zero 
calculation  is  skipped;  if  negative,  the  initial  density 
(Card  3.2),  is  taken  to  be  the  reference  density. 

10 

71-80 

THUG  -  Initial  Temperature  for  Hugoniot  Calculation.  If 
zero,  calculation  is  skipped;  if  negative,  the  initial 
temperature,  T^  (Card  3.2),  is  taken  to  be  the  reference 
temperature. 

Card  3.2 

-  FORMAT  (8E10.3) 

Variable 

Columns 

Description 

1 

1-10 

Number 

of  Elements  in  Material 

2 

11-20 

Switch  for  Type  of  Equation-of-State.  Must  be  -4  for 
AFWLE  equation-of-state. 

3 

21-30 

Pq  -  Reference  Density 

4 

31-40 

Tq  -  Reference  Temperature.  If  T^^  <  0,  code  sets  T^  = 
0.02567785  eV  (298®K). 

5 

41-50 

Blank 

6 

51-60 

Blank 

* 

7 

61-70 

r  -  Reference  Grtlneisen  Coefficient 

O 

8 

71-80 

Blank 

Card  3.3 

-  FORMAT  (8E10.3) 

Variable 

Columns 

Description 

9 

1-10 

Pmin  -  Lowest  Allowed  Solid  Density.  If  <  0,  code 

sets  p„,.„  =  0.8  p„. 
mm  0 

10 

11-20 

Blank 

11 

21-30 

Cy  -  Heat  Capacity.  If  C^  <  0,  code  sets  C^  =  3NoK. 

12 

31-40 

Eg  -  Sublimation  Energy  (cgs  units) 

13 

41-50 

V|^  -  Loading  Poisson's  Ratio 

14 

51-60 

V  -  Unloading  Poisson's  Ratio 

15 

61-70 

CL  -1 

Constants  Used  to  Describe  Density/Pressure 

16 

71-80 

CL  - 

“  Relationship  in  AFWLE 
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Card  3.4  -  FORMAT  (8E10.3) 


Variable 

Columns 

Description 

17 

1-10 

CU 

18 

11-20 

CZ 

19 

21-30 

P. 

_  Constants  Used  to  Describe  Density/ 

20 

31-40 

AM2 

Pressure  Relationship  in  AFWLE 

21 

41-50 

AM3 

22 

51-60 

BKM 

23  ■ 

61-70 

AMS  =  MUSTAR 

24 

71-80 

Blank 

Card  3.5 

-  FORMAT 

(5(F5.0,E10.3)) 

There  is 

one  pair 

of  the  following 

variables  for  each  element  in  variable 

of  card  3 

.2. 

1=1,  number  of  elements 

Odd  Variable  Z(I)  -  Atomic  Number  of  Element 

Even  Variable  Unnormalized  Atomic  Number  Fraction  of  Element  or  Unnormalized 
Atomic  Weight  Fraction  of  Element.  All  elements  should  be  de¬ 
fined  in  the  same  way. 


CARD  SET  12:  YIELD  CRITERIA 

When  the  elastic/plastic  option  is  to  be  used,  the  following  card  must  be  in¬ 
cluded  in  the  CSQ  II  input  card  stream.  There  must  be  one  card  for  each  AFWLE 
type  material.  When  the  material  is  to  be  treated  hydrodynamically,  a  blank 
card  must  be  entered. 

Card  12.1  -  FORMAT  (8E10.3) 


Variable 

Columns 

Description 

1 

1-10 

VM2  (dynes/cm^) 

2-4 

11-40 

Blank 

_ Constants  Used  to  Describe  AFWLE 

5 

41-50 

ST2 

"Yield  Curve 

6 

51-60 

Y2 

7 

61-70 

S2 

26 


Variable 

8 


f 


1 


Columns  Description 

71-80  AFWLE  Flag.  Must  be  0  to  treat  variables  as  AFWLE 
variables;  if  flag  =  0,  variables  will  receive  a 
CSQ  II  interpretation. 
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ABBREVIATIONS,  ACRONYMS.  AND  SYMBOLS 


C|_  elastic  wave  speed 

Cq  bulk  sound  speed 

Cy  heat  capacity  of  material 

E|^  mechtnlcal  energy  due  to  compression  (or  expansion)  on  Hugoniot 

H  Y  -  1 

K  Boltzmann's  constant 

N  Cj/re^ 

No  number  of  atoms  per  unit  mass 

P  pressure 

pressure  along  Hugoniot 
T  tracer  particle 

fj(p)  function  relating  pressure  to  density 

t  temperature 

t^  temperature  on  Hugoniot 

t^  ambient  temperature 

r  Grtfneisen  ratio  (material  constant) 

Y  ratio  of  specific  heats  of  the  perfect  gas 

e  internal  energy 

Cg  sublimation  energy  for  material 

n  1  -  Pq/p 

y  p/pQ  -  1 

V  Poisson's  ratio 

p  density 
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